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Abstract 



We have calculated the equation of state and elastic properties of face-centered cubic Fe and Fe-rich FcMg 
alloy at ultrahigh pressures from first principles using the Exact Muffin-Tin Orbitals method. The results 
show that adding Mg into Fe influences strongly the equation of state, and cause a large degree of softening 
of the elastic constants, even at concentrations as small as 1-2 at. %. Moreover, the elastic anisotropy 
increases, and the effect is higher at higher pressures. 
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Apart from a better understanding of one of our most 
important metals, one aim is to find the composition 
and structure of the Earth's core, which is likely to 
be made up of iron mixed with lighter elements. 

For a light element to be considered as a possible 
component in the Earth's core, it should naturally al- 
loy with Fe and be abundant within the Ea rth. While 
it is well-known that Mg is abundant ( Andersonl 
Il989h . until recently, magnesium was not considered 
as a possible candidate to be one of the light elements 
in the core, mainly because iron and magnesium differ 
too much in size to form alloys. However, it has been 
shown both theoretically and experimentally that at 
high pressures, the size mismatch of iron and mag- 
nesium h as decreased enough so that FeMg alloys 
can f orm ( Dubrovinskaia et al. , 2005: Tak afuii et al 



it is important to investigate elastic properties 
for the face-centered cubic structure (fee) of Fe 
rich FeMg a l loy. In particular, in a recent study 
( Asker et al. , 20091 ) we have demonstrated that the 
elastic anisotropy of fee Fe , estim ated according to 
Hill's definition (|Grimvalll . Il999h . is much higher 
than for hep Fe, and the effect persists in FeNi 
alloys. 

We report in this letter first-principles calculations 
of the equation of state and elastic properties of 
fee Fe, Fe 98 Mg Q2 , Fe 95 Mg 05 and Fe 9Q Mg 1Q at ul- 
trahigh pressures. The calculations have been car- 



20051 ) . At oxygen fugacity corresponding to Mg-MgO 



buffer, up to 10 at.% of Mg cou ld be dissolved in hep 



Fe at about 100 GPa pressure (jDubrovinskaia et al 



20051 ). In presence of magnesiowiistite (Mg,Fe)0 and 
dense silicate perovskite (M g, Fe)SiOa, about 1 a t.% 
Mg was observed in hep Fe (|Takafuii et alll2005[) . 

We acknowledge that the relatively high oxygen 
fugacity in the Earth interior, makes Mg a somewhat 
unlike ly candidate for the main light element in the 



core ( Dubrovinskaia et al. , l2005h . Still its presence 



has not been ruled out, and the study of its effect on 
the properties of Fe is motivated. 

The elastic properties of iron rich FeMg alloys were 
recentl y reported f o r the hexagonal close-packed 
(hep) (iKadas et all l2008h and for body-centered 
cubic (bec) (jKadas et all [2009) structures. How- 
ever, Mikhaylushkin et al. have shown that the 
fee phase of Fe can not be ruled out as a possible 
structure stable at conditions corresp o nding to the 



ried o ut using the EMTO-CPA method (jVitos et al 
200l[ ) (described below), which is a suitable tool for 
calculating elastic properties in alloys. 

The paper is organized as follows: in Sec.[2]we out- 
line the theory used for calculating elastic properties 
from first-principles by the EMTO-CPA method to- 
gether with details of our calculations and in Sec. 
we present the equation of state, the single crystal 
and polycrystalline elastic properties of fee FeMg. 

2 Theory 

In this work, we first calculated the total energy for a 
set of atomic volumes between 41.6, and 78.0 Bohr 3 . 
The results from these calculations were then used 
to find the equation of s tate (EOS) through the 
Birch-Murnaghan scheme (jBirchl . Il947l ) (3rd order). 
The EOS provides a relation between volume and 
pressure as well as allows us to study the pressure 
dependence of the bulk modulus. 

In cubic structures, there are th ree independent 
elastic constants: en, C12 and C44. ( Grimvall . 19991 ) 
The c' elastic constant (also known as Zener's elas- 
tic constant) can be obtained from the orthorhom- 
bic st ructure with the following deformation I Vitosl . 
20071 ): 



Earth's core. (jMikhavlushkin et all [2007) Since 
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The change in total energy upon distortion is 

AE = 2V ■ d ■ S 2 + 0(8 4 ) (2) 

Now the en and c\i can be calculated from the 
relations between c! and the bulk modulus: 
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From Eqs. (|2I3I4[) the c\\ and C\i elastic constants 
can be obtained. 

Next, the C44 elastic constant c an be obtain ed by 
performing monoclinic distortion ( Vitod . 2007 ): 
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In this case the change in energy is: 

AE = 2V • c 44 • S 2 + £>(<5 4 ) (6) 

In the calculations of both c' and C44 the dis- 
tortions used were 5 = [0.00, 0.01, 0.05]. The 
elastic constants were then obtained from (linear) 
interpolations of the total energy as function of the 
square of the distortion, AE tot = AE tot (S 2 ). The 
motivation for this procedure is that it simplifies 
the choice of suitable distortions, and is numerically 
preferable over interpolation of distortion. In this 
work we only do volume-conserving distortions. 
The reason for this is that we seek to calculate 
the elastic properties as functions of pressure, but 
in the computational setup we must use volume 
(rather than pressure) as input parameter. The EOS 
can then be used to obtain the elastic constants 
as functions of pressure. Also, the total energy 
depends more strongly on volume than on the dis- 
tortions ( Steinle-Neumann et al. . 19991 ). indicating 
that it should be more stable numerically to do 
volume-conserving distortions whenever possible. 

The elastic constants described above (except B) 
are single-crystal elastic constants. When doing mea- 
surements, one is often dealing with polycrystalline 



samples of a material, and hence the single-crystal 
constants should be averaged in a suitable way to 
facilitate comparison between experiments and theo- 
retical calculations. The polycry stalline shear mo duli 
defined by Reuss and Voigt are (lGrimval]| . [T999h : 
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respectively. In the case of cubic systems, the Reuss 
and Voigt bulk moduli are equal, and defined in 
Eq.©. 

Further, the elastic anisotropy according to Hill is 
the weighted difference between the Voigt and Reuss 
shear moduli: 



wrh — 



Gy — Gr 
Gy + Gr 



(9) 



and hence is called the Voigt- Reuss -Hill (VRH) 
anisotropy (|Grimvalj[l999l) . 



All the calculations in this work have been 
carrie d out within Density-Functional The- 
ory (jHohenberg and Kohnl 119641: iKohn and Shaml . 
19651 ) (DFT) using the Exact Muffin-t in Or- 
(EMTO) 



( Andersen et all . 11994 : 



bital Metho d 

Vitosl . Il999l l200lh . The problem of disorder was 
treated with the Coherent Po tential App roxima- 



tion (CPA) l|Sovenl . 119671 : IVelickv et all 1968 



implemented in the EMTO-CPA (jVitos et aU l200ll 
method. For exchange and correlation the Local- 
Density approximation (LDA) , parametrized b y 



1996). 



Perdew, Burke and Ernzerhof (jPerdew et al. 
was used. After self-consistency had been reached 
in the LDA calculations, the total energy was calcu- 
lated with Full Charge Density Technique (|Vitosl . 



2007; Vitos et al 



19971) (FCD). In this last step the 
Generalized Gradient Approximation (GGA) was 
used. This procedure provides essentially the same 
materials properties as those of a self-consistent 
GGA scheme for non-magnetic systems. While 
calculations are done for T = K, the effect of 
temperature on elastic properti es is expected to 
decrease with increasing pressure ( Anders on. 1989). 



2 



We have performed the calculations for Fe, 
Fe ss Mg 02 , Fe 95 Mg 05 and Fe 90 Mg w in order to 
investigate the effects of alloying Mg in Fe at small 
Mg content. 

For the EMTO calculations, the integration in re- 
ciprocal space was performed over a grid of 29x29x29 
points, and the Green's function was integrated using 
16 points in the complex energy plane. The basis set 
used s,p,d,f orbitals. 

3 Results 

The volumes as function of pressure for pure Fe and 
for FegoMgio are shown in Fig. [TJ From this figure 
it is clear that adding Mg to Fe increases the volume 
corresponding to a certain pressure. This is expected 
due to the larger atomic size of Mg, and is well-known 
from earlier work. When the pressure is increased, 
the difference in volumes decreases, which is due to 
the fact that Mg has higher compressibility than Fe. 
( Dubrovinskaia et all [2005) 



loys (|Kadas et all |2008( ) . The density for Fe go Mg 10 
lies close to those of the Preliminary Referenc e Eart h 
Model (PREM) (jPziewonski and Andersonl . Il98ll ). 
However, including temperature effects in our cal- 
culations is likely to change the agreement. From 
the inset, showing the concentration dependence at 
P = 350 GPa, we see that adding 2 % Mg causes the 
density to decrease about 1.5 %. 
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Figure 1: (Color online) Volume as function of pres- 
sure in fee Fe^Mgio (red full line) and pure Fe (black 
dashed line). 

In Fig. [5] we report the density as function of pres- 
sure. The results show that adding Mg into Fe de- 
creases the density, which is expected. The results 
are similar to that reported earlier for hep Fe-Mg al- 



Figure 2: (Color online) Density as function of pres- 
sure in fee FegoMgio (red full line) and pure Fe (black 
dashed line) . The inset shows the density as function 
of Mg content for P = 350 GPa. 

Next, Fig. [3] shows the bulk modulus as function of 
pressure. Also here the results are sim ilar to those 
reported for hep (jKadas et all l2008h . The bulk 



modulus is about 55 GPa lower for FegoMg w than 
for pure Fe over the whole pressure range consid- 
ered here, indicating that the FeMg alloy is indeed 
softer than pure Fe. Also here we note that the 
bulk modu lus for FegoMgio is c lose to that of the 



PREM (jPziewonski and Anderson! . 119811 ). but tem- 
perature effects may change this, and other elastic 
properties may not agree to such a large degree. Still, 
this work aims to show the effects of adding Mg in Fe, 
which is still likely to be important at elevated tem- 
peratures. At 350 GPa, adding 2 % Mg decreases the 
bulk modulus by about 1.0 %. 

The single- crystal elastic constants are shown in 
Fig. [H It is interesting to note that both en and C44 
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decrease Gy by 4.3 % and Gr by 6.9 %. 

Fig.[6]shows the Voigt-Reuss-Hill elastic anisotropy 
(Avrh) as function of pressure for pure fee Fe and 
for fee FegoMgio- This graph shows that adding 
10% Mg into Fe cause a remarkable increase of the 
anisotropy. At 300 GPa, the anisotropy is about 40% 
higher for fee FegoMgio than for pure fee Fe. The 
pressure dependence of the anisotropy is also higher 
for FegoMgio than for pure Fe. 



4 Conclusion 



Figure 3: (Color online) Bulk modulus as function 
of pressure in fee FegoMgio (red full line) and pure 
Fe (black dashed line). The inset shows the bulk 
modulus as function of Mg content for P = 350 GPa. 

show considerable softening when adding Mg. This 
is also evident for c', which is shown for comparison. 
However, the cyi elastic constant is almost constant, 
and for very high pressures it even increases a little 
when adding Mg. This can be understood from the 
relations used to obtain en and C12 from c' and the 
bulk modulus: 

4 

en =B+-c' 

C12 =B- y (io) 

Since both the bulk modulus and c' show similar be- 
havior when adding Mg, en will follow the same 
trend. However, Eq. (fT0|) suggests that for C12 the 
effect is to some extent canceled out. This can also 
be seen from the inset for en in Fig. 0] For P = 350 
GPa, adding 2 % Mg decreases c' by about 8.2 %, en 
by about 2.4 %, ci 2 only 0.1 % and C44 by 3.6 %. 

The polycrystallinc clastic constants are shown in 
Fig. This figure shows that not just the single- 
crystal elastic constants change considerably when 
going from pure Fe to FegoMg w . A large degree 
of softening takes place for both Gr and Gy when 
adding Mg to Fe. For P = 350 GPa, adding 2 % Mg 



We have calculated the elastic properties of fee Fe 
and FegoMgio random alloy up to pressures of the 
Earth's core from first-principles methods. We show 
that adding 10 at. % Mg into iron has a substantial 
effect on the single crystal elastic constants. More- 
over, even as small Mg concentrations as 2 % sub- 
stantially affect some elastic constants, ie c' which 
decreases by 8.2 % at 350 GPa. Further we find that 
the effect on the elastic anisotropy is remarkably high 
and is higher for the fee FeMg alloy than for pure fee 
Fe. We therefore conclude that the possibility of Mg 
in the Earth's core must be included in models, even 
though it may not be the dominant light element. 

5 Acknowledgments 

Useful discussions with L. Vitos are gratefully ac- 
knowledged. We are grateful for financial support 
from the Swedish research council and the Goran 
Gustafsson Foundation for Research in Natural Sci- 
ences and Medicine. The National Supercomputer 
centre (NSC) and Center for Parallel Computers 
(PDC) are acknowledged for computer support. 

References 

Andersen, O. K., Jepsen, O., Krier, G., 1994. Exact 
Muffin-Tin Orbital Theory. World Science Publish- 
ing Co. 



4 



Anderson, D. L., 1989. Theory of the Earth. Black- 
well Scientific Publications. 

Asker, C, Vitos, L., Abrikosov, I. A., 2009. Elas- 
tic constants and anisotropy in feni alloys at high 
pressures from first-principles calculations. Phys- 
ical Review B (Condensed Matter and Materials 
Physics) 79 (21), 214112. 

Bclonoshko, A. B., Ahuja, R., Johansson, B., August 
2003. Stability of the body-centered cubic phase of 
iron in the Earth's inner core. Nature 424, 1032. 

Bclonoshko, A. B., Skorodumova, N. V., Rosengren, 
A., Johansson, B., 2008. Elastic Anisotropy of 
Earth's Inner Core. Science 319 (5864), 797-800. 

Birch, F., Jun 1947. Finite clastic strain of cubic crys- 
tals. Phys. Rev. 71 (11), 809-824. 

Dubrovinskaia, N., Dubrovinsky, L., Kantor, I., 
Crichton, W. A., Dmitriev, V., Prakapenka, V., 
Shcn, G., Vitos, L., Ahuja, R., Johansson, B., 
Abrikosov, I. A., 2005. Beating the miscibility bar- 
rier between iron group elements and magnesium 
by high-pressure alloying. Physical Review Letters 
95 (24), 245502. 

Dubrovinsky, L., Dubrovinskaia, N., Langenhorst, F., 
Dobson, D., Rubie, D., Geszmann, C, Abrikosov, 
I. A., Johansson, B., Baykov, V. I., Vitos, L., Bi- 
han, T. L., Crichton, W. A., Dmitriev, V., Weber, 
H.-P., Mar. 2003. Iron-silica interaction at extreme 
conditions and the electrically conducting layer at 
the base of Earth's mantle. Nature 422 (6927), 58- 
61. 

Dubrovinsky, L., Dubrovinskaia, N., Narygina, O., 
Kantor, I., Kuznetzov, A., Prakapenka, V. B., 
Vitos, L., Johansson, B., Mikhaylushkin, A. S., 
Simak, S. I., Abrikosov, I. A., Jun. 2007. Body- 
Centered Cubic Iron-Nickel Alloy in Earth's Core. 
Science 316 (5833), 1880-1883. 

Dzicwonski, A. M., Anderson, D. L., Jun. 1981. Pre- 
liminary reference earth model. Physics of The 
Earth and Planetary Interiors 25 (4), 297-356. 



Grimvall, C, 1999. Thcrmophysical Properties of 
Materials, 1st Edition. Elsevier. 

Hohenberg, P., Kohn, W., November 1964. Inhomo- 
geneous Electron Gas. Phys. Rev. 136 (3B), B864- 
B871. 

Kadas, K., Vitos, L., Ahuja, R., Jul. 2008. Elas- 
tic properties of iron-rich hep Fe-Mg alloys up to 
Earth's core pressures. Earth and Planetary Sci- 
ence Letters 271 (1-4), 221-225. 

Kadas, K., Vitos, L., Johansson, B., Ahuja, R., 2009. 
Stability of body-centered cubic iron-magnesium 
alloys in the Earth's inner core. Proceedings of 
the National Academy of Sciences 106 (37), 15560- 
15562. 

Kohn, W., Sham, L. J., November 1965. Self- 
Consistent Equations Including Exchange and Cor- 
relation Effects. Phys. Rev. 140 (4A), A1133- 
A1138. 

Kuwayama, Y., Hirose, K., Sata, N., Ohishi, Y., 2008. 
Phase relations of iron and iron-nickel alloys up to 
300 GPa: Implications for composition and struc- 
ture of the Earth's inner core. Earth and Planetary 
Science Letters 273 (3-4), 379 - 385. 

Lin, J.-F., Heinz, D. L., Campbell, A. J., Devine, 
J. M., Shen, G., Jan. 2002. Iron-Silicon Alloy in 
Earth's Core? Science 295 (5553), 313-315. 

Mikhaylushkin, A. S., Simak, S. I., Dubrovinsky, L., 
Dubrovinskaia, N., Johansson, B., Abrikosov, I. A., 
Oct. 2007. Pure iron compressed and heated to ex- 
treme conditions. Phys. Rev. Lett. 99 (16), 165505- 
4. 

Perdcw, J. P., Burke, K., Ernzerhof, M., Oct 1996. 
Generalized gradient approximation made simple. 
Phys. Rev. Lett. 77 (18), 3865-3868. 

Soven, P., April 1967. Coherent-Potential Model 
of Substitutional Disordered Alloys. Phys. Rev. 
156 (3), 809-813. 

Steinle-Neumann, G., Stixrude, L., Cohen, R. E., Jul 
1999. First-principles elastic constants for the hep 



5 



transition metals Fe, Co, and Re at high pressure. 
Phys. Rev. B 60 (2), 791-799. 

Steinle-Neumann, G., Stixrude, L., Cohen, R. E., 
Gulscren, O., Sep. 2001. Elasticity of iron at the 
temperature of the Earth's inner core. Nature 
413 (6851), 57-60. 

Takafuji, N., Hirose, K., Mitome, M., Bando, Y., 
Mar. 2005. Solubilities of o and si in liquid iron 
in equilibrium with (mg,fc)sio3 perovskite and the 
light elements in the core. Geophys. Res. Lett. 32, 
L06313. 

Vclicky, B., Kirkpatrick, S., Ehrenreich, H., 1968. 
Second- and Third- order Elastic Constants of Al- 
kali Metals. Phys Rev. 147, 747. 

Vitos, L., November 1999. Application of the exact 
muffin-tin orbital theory: the sperical cell approx- 
imation. Comp. Mat. Sci. 18, 24. 

Vitos, L., Jun 2001. Total-energy method based on 
the exact muffin-tin orbitals theory. Phys. Rev. B 
64 (1), 014107. 

Vitos, L., 2007. Computational Quantum Mechanics 
for Materials Science The EMTO Method and Ap- 
plications, 1st Edition. Springer Verlag. 

Vitos, L., Abrikosov, I. A., Johansson, B., Sep 2001. 
Anisotropic Lattice Distortions in Random Al- 
loys from First-Principles Theory. Phys. Rev. Lett. 
87 (15), 156401. 

Vitos, L., Kollar, J., Skriver, H. L., May 1997. Full 
charge-density scheme with a kinetic-energy cor- 
rection: Application to ground-state properties of 
the 4d metals. Phys. Rev. B 55 (20), 13521-13527. 

Vocadlo, L., 2007. Ab initio calculations of iron and 
iron alloys at inner core conditions: Evidence for a 
partially molten core? Earth Planet Sci Lett 254, 
227. 

Vocadlo, L., Alfe, D., Gillan, M. J., Wood, I. G., 
Brodholt, J. P., Price, G. D., July 2003. Nature 
424, 536. 



Vocadlo, L., Wood, I. G., Alfe, D., Price, G., Apr. 
2008. Ab initio calculations on the free energy 
and high P-T elasticity of face-centred-cubic iron. 
Earth and Planetary Science Letters 268 (3-4), 
444-449. 



6 




100 150 200 250 300 350 
Pressure (GPa) 

Figure 4: (Color online) Elastic constants as function 
of pressure in fee FegoMg w (red full line) and pure Fe 
(black dashed line) . All elastic constants are in GPa. 
The inset in each subfigure show the corresponding 
elastic constant as function of Mg content for P = 
350 GPa. Note that all the insets have the same scale 
(250 GPa). 
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Figure 5: (Color online) Polycrystalline shear modu- 
lus as function of pressure. Black lines and squares 
denote results for fee Fe, while red lines and circles 
denote fee FegoMgiQ. Further, full lines and filled 
symbols denote Gr, while dashed lines and empty 
symbols denote Gy. The inset shows Gy (open dia- 
monds) and Gr (filled diamonds) as function of Mg 
content for P = 350 GPa. 



I Fe 

> Fe 9 o M 9i 




Mg cone (at. %) 



100 150 200 
Pressure (GPa) 



Figure 6: (Color online) Voigt-Reuss-Hill anisotropy 
as function of pressure in fee FegoMgio (red full line) 
and pure Fe (black dashed line). The inset shows 
the anisotropy as function of Mg content for P = 
350 GPa. 
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